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Abstract 

The intimate geometry of space-time is expected to suffer stochastic fluctuations as 
a result of quantum gravitational effects. These fluctuations may induce observable 
consequences on the propagation of high energy particles over large distances, so that 
the strength and the characteristics of these fluctuations may be constrained, mainly 
in the range of energies of interest for cosmic ray physics. While invoked as a possible 
explanation for the detection of the puzzling cosmic rays with energies in excess 
of the threshold for photopion production (the so-called super-GZK particles), we 
demonstrate here that lower energy observations may provide strong constraints on 
the role of a fluctuating space-time structure. 



1 Introduction 



In the recent past it has been pointed out [1-4,6] that observations in cosmic 
ray physics, and in particular Ultra High Energy Cosmic Ray (UHECR) ex- 
periments can be used to constrain to high precision the possibility that the 
fundamental Lorentz invariance (LI) may be broken to some level. 

In fact, violations of LI have been invoked by many authors [3-5] to explain 
puzzling observations related to the detection of cosmic rays with energy in 
excess of ~ 10 20 eV, the so-called Greisen-Zatsepin-Kuzmin (GZK) [7] cut- 
off, and the more uncertain excess of photons in the TeV energy region from 
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distant sources. The two phenomena have in common that in both cases the 
detected particles should instead suffer significant attenuation since their en- 
ergy is above the threshold for photopion production (in the first case) and 
pair production [8] (in the second case). In particular the GZK cutoff has been 
expected for a long time because of the dramatic energy losses of high energy 
protons with the 2>°K Cosmic Microwave Background Radiation (CMBR). Vi- 
olations of LI in both cases may act on the kinematics of the processes so to 
move the thresholds in regions which are possibly outside the observational 
window. In the following we mainly concentrate upon the case of UHECRs, 
pointing out, when necessary, the differences with respect to the case of pair 
production. 

The recipes for the violations of LI that have been suggested insofar generally 
consist of requiring an explicit modification of the dispersion relation of UHE 
particles, due to their propagation in the "vacuum" , now affected by quantum 
gravity (QG). This effect is generally parametrized by introducing a typical 
mass, expected to be of the order of the Plank mass (Mp), that sets the scale 
for QG to become effective. 

In fact, explicit modifications of the dispersion relation are not really neces- 
sary, as was recently pointed out in Refs. [9-12] for the case of propagation 
of UHECRs. There is an intrinsic uncertainty in the process of measurement 
of physical quantities on time (space) scales comparable with the Planck time 
(length). It is generally believed that coordinate measurements cannot be per- 
formed with precision better than Planck distance (time) 5x > lp, since such a 
measurement would result in the production of a black hole. This means that 
the metric of space-time must feature quantum fluctuations on the Planck 
scale. A similar line of thought implies that an uncertainty in the measure- 
ment of energy and momentum of particles must occur, as described by the 
relation 5p ~ 8E ~ E 2 /M P . According to the authors of Refs. [10,11] the ap- 
parent problem of super-GZK particles may find a solution also in the context 
of this uncertainty approach. 

We discuss here this appealing approach more in detail, by taking into account 
the effects of the propagation of CRs in the QG vacuum in the presence of 
the universal microwave background radiation. A fluctuating metric implies 
that different measurements of the particle energy or momentum may result 
in different outcomes. Therefore it becomes important to define the probabil- 
ity that the measured energy (momentum) of a particle is above some fixed 
value. Note that averaging over a large number of measurements would yield 
the classical values for the energy and momentum. The process of measure- 
ment mentioned above, during the propagation of particles over cosmological 
distances occurs at each single interaction of the particle with the environ- 
ment. At each interaction vertex, the fluctuating energy /momentum of the 
particle is compared with the kinematic threshold for the occurrence of some 
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physical process (in our case the photopion production). A clear consequence 
of this approach is that particles with classical energy below the standard 
Lorentz invariant threshold have a certain probability of interacting. In the 
same way, particles above the classical threshold have a finite probability of 
evading interaction. 

We show here that the most striking consequences of the approach described 
above derive from low energy particles rather than from particles otherwise 
above the threshold for photopion production. 



2 The effect of Space-Time fluctuations on propagation of UHE 
particles. 

While electroweak and strong interactions propagate through space-time, grav- 
ity turns out to be a property of the space-time itself. This simple statement 
has profound implications in that our belief that gravity can be turned into 
a quantum theory immediately implies that the structure of space-time has 
quantum fluctuations. Another way of rephrasing this concept is that space- 
time is expected to have a granular (or foamy) structure, where however the 
size of space-time cells fluctuates stochastically, thereby causing an intrinsic 
uncertainty in the measurements of space-time lengths, and indirectly of en- 
ergy and momentum of a particle moving through space-time. The uncertainty 
appears on scales comparable with the Planck scale. 

It is generally argued that measurements of distances (times) smaller than 
the Planck length (time) are conceptually unfeasible, since the process of 
measurement collects in a Planck size cell an energy in excess of the Planck 
mass, hence forming a black hole, in which information is lost. This can be 
traslated in different ways into an uncertainty on energy-momentum measure- 
ments ([10,11]). The Planck length is a good estimate of the uncertainty in the 
De Broglie wave-length A of a particle with momentum p. Therefore 5 A ~ lp, 
and 5p = 5(1/ A) « (p 2 l P ) = (p 2 /M P ). 

Speculating on the exact characteristics of the fluctuations induced by QG 
is beyond the scope of the present paper, and it would probably be useless 
anyway, since the current status of QG approaches does not allow such a kind 
of knowledge. We decided then to adopt a purely phenomenological approach, 
in which some reasonable assumptions are made concerning the fluctuations 
in the fabric of space-time, and their consequences for the propagation of 
high energy particles are inferred. Comparison with experimental data then 
possibly constrains QG models. 

Following [10], we assume that in each measurement: 
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the values of energy (momentum) fluctuate around their average values 
(assumed to be the result theoretically recoverable for an infinite number of 
measurements of the same observable): 

E 2 

E^ E + a 



M P 

2 



P"P + />^ (2) 

with a, (3 normally distributed variables and p the modulus of the 3-momentum 
(for simplicity we assume rotationally invariant fluctuations); 
• the dispersion relation fluctuates as follows: 

P^P u = E 2 -p 2 + 1 j^- = m 2 (3) 

and 7 is again a normally distributed variable. 

Ideally, QG should predict the type of fluctuations introduced above, but, as 
already stressed, this is currently out of reach, therefore we assume here that 
the fluctuations are gaussian. Our conclusions are however not sensitive to 
this assumption: essentially any symmetrical distribution with variance « 1, 
within a large factor, would give essentially the same results. Furthermore 
we assume that a, (3 and 7 are uncorrelated random variables; again, this 
assumption reflects our ignorance in the dynamics of QG Q. 

Our interest will be now concentrated upon processes of the type 



a + b — ► c + d 



where we assume that a kinematic threshold is present; in the realm of UHECR 
physics (a,b) is either (7,73^) or (p, j^k) an d (c,d) is (e + , e~) or (N,tt). 

To find the value of initial momenta for which the reaction occurs we write 
down energy-momentum conservation equations and solve them with the help 
of the dispersion relations, as discussed in detail in [6]. 

The energy momentum conservation relations are (in the laboratory frame, 
and specializing to the case in which the target (b) is a low energy background 



The fluctuations described above will in general derive from metric fluctuations 
of magnitude 5g^ u ~ h^ u ^j- [3,11]. Our assumption reflects the fact that, while the 
magnitude of the fluctuation can be guessed, we do not make any assumption on 
its tensorial structure W . 
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photon for which fluctuations can be entirely neglected) 



These equations refer to head-on collisions and collinear reaction products, 
which is appropriate for threshold computations. Together with the modified 
dispersion relations, these equations, after some manipulations, lead to a cubic 
equation for the initial momentum as a function of the momentum of one of 
products, and, after minimization, they define the threshold for the process 
considered. In figure 1 we report the distribution of thresholds in the ~ 70% 
of cases in which the solution is physical; in the other cases the kinematics 
does not allow the reaction. 




10 

Energy (e^ 



Fig. 1. Threshold distribution for p^z°K Ntt. In the 30% of cases the reaction is 
not allowed. 

This threshold distribution can be interpreted in the following way: a particle 
with energy above ~ 10 15 eV has essentially 70% probability of being above 
threshold, and therefore to be absorbed. In the other 30% of the cases the 
protons do not interact. 

In (4,5) the fluctuations are taken independently for each particle, which is 
justified as long as the energies are appreciably smaller than the Planck energy. 
At that point it becomes plausible that different particles experience the same 
fluctuations, or more precisely fluctuations of the same region of space-time. 
It is instructive to consider this case in some more detail: we introduce then 
the four-momenta (and dispersion relations) of all particles fluctuating in the 
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same way. Specializing to proton interaction on CMBR, the equation which 
defines the threshold p t h is [6] : 



2pg m^rrip 

(ml + 2m 7T m p )M P (m n + m p )' A 




(6) 



where 77 is a gaussian variable with zero average and variance of the order 
of (but not exactly equal to) one, and p is the L.I. threshold (GZK). The 
threshold is the positive solution of this equation. 

The coefficient of the cubic term is very large, of the order of 10 13 in this case, 
so that unless 77 is 0(1CT 13 ), we can write, neglecting pion mass 

f m 2 p M P y 

Pth ~ Po 3 • 7 
V VPo J 



When 77 becomes negative, the above equation has no positive root; this hap- 
pens essentially in 50% of the cases. Since the gaussian distribution is flat 
in a small interval around zero, the distribution of thresholds for positive 77 
peaks around the value for 77 1, meaning that the threshold moves almost 
always down to a value of « 10 15 eV [6]; essentially the same result holds 
for fluctuations affecting only the incident (higest energy) particle. For inde- 
pendent fluctuations of final momenta, the asymmetry in the probability of 
interaction arises from the fact that even exceedingly small negative values of 
the fluctuations lead to unphysical solutions. 

Building upon our findings, we now apply the same calculations to the case of 
UHECR protons propagating on cosmological distances. An additional ingre- 
dient is needed to complete the dynamics of the process of photopion produc- 
tion, namely the cross section. The rather strong assumption adopted here is 
that the cross section remains the same as the Lorentz invariant one, provided 
the reaction is kinematically allowed. This implies that the interaction lengths 
remain unchanged. 

In order to assess the situation of UHECRs, we first consider the case of 
particles above the threshold for photopion production in a Lorentz invariant 
world. According with eqs. (4, 5), in this case particles have a probability 
of ~ 30% of being not kinematically allowed to interact inelastically with 
a photon in the CMBR. Therefore, if our assumption on the invariance of 
the interaction length is correct, then each proton is still expected to make 
photopion production, although with a slightly larger pathlength. 

The situation is however even more interesting for particles that are below 
the Lorentz invariant threshold for the process of photopion production. If the 
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energy is below a few 10 18 eV, a galactic origin seems to be in good agreement 
with measurements of the anisotropy of cosmic ray arrival directions [13,14]. 
We will not consider these energies any longer. On the other hand, at energies 
in excess of 10 19 eV, cosmic rays are believed to be extragalactic protons, 
mainly on the ground of the comparison of the size of the magnetized region 
of our Galaxy and the Larmor radius of these particles. We take these pieces 
of information as the basis for our line of thought. If the cosmic rays observed 
in the energy range E > 10 19 eV are extragalactic protons, then our previous 
calculations apply and we may expect that these particles have a ~ 70% 
probability of suffering photopion production, even if their energy is below the 
classical threshold for this process. Note that the pathlength associated with 
the process is of the order of the typical pathlength for photopion production 
(a few tens of Mpc), therefore we are here discussing a dramatic process in 
which the absorption length of particles drops from Gpc, which would be 
pertinent to particles with energy below ~ 10 20 eV in a Lorentz invariant 
world, to several Mpc, with a corresponding suppression of the flux. What 
are the consequences for the observed fluxes of cosmic rays? The above result 
implies that all with E > 10 15 eV are produced within a radius of several 
tens of Mpc, and above this energy there is no dramatic change of pathlength 
with energy. There is no longer anything special about E ~ 10 20 eV, and any 
mechanism invoked to explain the flux of super-GZK particles must be at work 
also at lower energies. 

The basic situation remains the same in the case of pair production as the 
physical process under consideration. For a source at cosmological distance, 
a cutoff is expected due to pair production off the far infrared background 
(FIR) or the microwave background. Using the results in [6] we expect that 
the modified thresholds are a factor 0.06 (0.73) lower than the Lorentz in- 
variant ones for the case of interaction on the CMBR (FIR). There is also a 
small increase in the pathlengths above the threshold, which would appear 
exponentially in the expression for the flux. Therefore there are two effects 
that go in opposite directions: the first moves the threshold to even lower en- 
ergies, and the second increases the flux of radiation at Earth because of the 
increase of the pathlength. It seems that geometry fluctuations do not provide 
an immediate explanation of the possible detection of particles in excess of 
the expected ones from distance sources in the TeV region. In any case the 
experimental evidence for such an excess seems at present all but established. 



3 Discussion, conclusions and perspectives. 

The investigation of the dependence of the kinematic thresholds for physical 
processes has been shown to be a powerful tool to study the Physics of space- 
time at extremely high energies, close to that Planck scale where the fabric of 
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space-time is expected to change, from the flat and well-behaved sheet that 
we experience in everyday life to a complicated and unpredictable foam of 
probability that, till now, evaded any kind of direct investigation. 

In this paper we related this foamy structure of space-time with the process of 
measurement of the physical properties of particles, in particular their energy 
and momentum. We found that the fluctuating metric may induce a violation 
of Lorentz invariance that changes the thresholds for the photopion production 
of a very high energy proton off the photons of the CMBR, or for the pair 
production of a high energy gamma ray in the bath of the FIR or CMBR 
photons. 

For the case of UHECRs interacting with the CMBR, we obtained a picture 
that changes radically our view of the effect of QG on this phenomenon, as 
introduced in previous papers: not only particles with energy above ~ 10 20 
eV are affected by the fluctuations in space-time, but also particles with lower 
energy, down to ~ 10 15 eV seem to be affected by such fluctuations. In fact 
the latter, as a result of a fluctuating space-time, may end up being above the 
threshold for photopion production, so that particles may suffer significant 
absorption. Our conclusion is that all particles with energy in excess of ~ 10 15 
eV eventually detected at Earth would be generated at distances comparable 
with the pathlength for photopion production (~ 100 Mpc). A consequence 
of this is that there is no longer anything special characterizing the energy 
~ 10 20 eV. 

Since the conclusion of our work are quite strong, it is important to summarize 
in detail the assumptions involved in the calculations: 

1) space-time fluctuations follow from Quantum Gravity. This is a rather mild 
statement and essentially accepted on the basis of the "space-time foam" ap- 
proach [15]. Following [10,11] we further assume that these fluctuations are at 
the level of the Planck scale 

2) Fluctuations in the momentum and in the dispersion relation are induced 
by the fluctuations in the metric of space-time. Our assumption is that the 
fluctuations in the dispersion relation and in the energy and momenta of the 
particles involved can be considered as independent. This is an assumption 
usually shared by most literature in the field. 

3) The cross section for photopion production at threshold is assumed to be the 
same as in the Lorentz invariant case. The cross section can be considered as 
a combination of a matrix element and a phase space factor. While the former 

2 We assume that this is true in the comoving reference frame, in which the CMB 
radiation is isotropic. If this is assumed to be true in all reference frames, than it is 
necessary to modify the Lorentz transformations [16,17] 
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may well be affected by violation of Lorentz invariance, it can be demonstrated 
that the phase space does not change appreciably once the incident particle 
has energy above the threshold for photopion production. This may be not 
true in other situations, as discussed in the literature [18,19]. 

We list below some tests that may allow to understand whether the current or 
future observations are compatible with the scenario discussed in this paper. 

a) Future experiments [20] dedicated to the detection of UHECRs will pro- 
vide a substantial increase in the statistics, so that the spectral features of the 
UHECRs in the energy region E > 10 19 eV can be resolved, and further indi- 
cations on the nature of primaries and their possible extragalactic origin will 
be obtained. In particular the present possible disgreement between AGASA 
[21] and HiRes [22] will be clarified. 

One should also keep in mind that an evaluation of the expected flux in terms 
of sources distributed as normal galaxies is in contradiction with AGASA 
data by an amount ranging from 2 to 6a depending on the assumed source 
spectrum [23]. Since the nature of the sources is not known, it is not clear if 
their abundance within the absorption pathlength is sufficient to explain the 
observed flux in presence of space-time fluctuations, nor if they can induce 
observable anisotropies. 

In any case, in a Lorentz invariant framework a suppression in the flux at 
~ 10 20 eV is expected. If such a feature is unambiguously detected in the 
UHECR spectrum, no much room would be left for the fluctuations of space- 
time discussed in this paper, since in this scenario nothing special happens 
around 10 20 eV. In quantitative terms [6] this would imply a phenomenological 
bound on l P now interpreted as a parameter: l p < 10~ 46 cm instead of l p ~ 
10 -33 cm; in other words, only fluctuations with variance ~ 10~ 13 , instead of 
1, would be allowed 

b) According with our findings, all particles with energy in excess of ~ 10 15 
eV lose their energy by photopion production on cosmological spatial scales, 
as a result of the metric fluctuations. This energy ends up mainly in gamma 
rays, neutrinos and protons. The protons pile up in the energy region right 
below ~ 10 15 eV. The gamma ray component actually generates an electro- 
magnetic cascade that ends up contributing low energy gamma rays, in the 
energy band accessible to instruments like EGRET and GLAST. This cascade 
flux cannot be larger than the measured electromagnetic energy density in the 
same band u^f = 10~ 6 eV/cm 3 [24]. The cascade flux in our scenario can be 
estimated as follows. Let $(£7) = &q(E/Eq) ~ 7 be the emissivity in UHECRs 

3 Alternatively, one can assume a more general form of fluctuations, i.e. 5E ~ 
E(E/Mp) a and similar for momentum and dispersion relations. In this case a bound 
a > 2.3 would follow. 
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(particles/cm 3 /s/GeV). Let us choose the energy E = 10 10 GeV and let us 
normalize the flux to the observations at the energy E . The total energy going 
into the cascade can be shown to be 

5 x 1(T 4 2 _„ , _ 3 

^cas ~ 7T X min t eV CUl , 

7-2 

where £ is the fraction of energy going into gamma rays in each photopion 
production, and x min = (E th /E ) = 10~ 4 for E th = 10 15 eV. It is easy to see 
that, for 7 = 2.7, the cascade bound is violated unless £ <C 1CT 3 . 

One note of warning has to be sent concerning the development of the elec- 
tromagnetic cascade: the same violations of LI discussed here affect other 
processes, as stressed in the paper. For instance pair production and pion de- 
cay are also affected by violations of LI [18]. Therefore the possibility that the 
cascade limit is exceeded concerns only those scenarios of violations of LI that 
do not inhibit appreciably pair production and the decay of neutral pions. 

The protons piled up at energies right below 10 15 eV, would be a nice signature 
of this scenario, but it seems difficult to envision a way of detecting these 
remnants. In fact, even a tiny magnetic field on cosmo logical scales would make 
the arrival time of these particles to Earth larger than the age of the universe. 
Moreover, even assuming an exactly zero extragalactic magnetic field, these 
particles need to penetrate the magnetic field of our own Galaxy and mix with 
the galactic cosmic rays, making their detection extremely problematic if not 
impossible. 

We believe that the results presented here are robust, in the sense that the 
qualitative features are rather insensitive to modifications of the theoretical 
scheme. Clearly a more detailed flux computation, taking into account prop- 
agation of primaries as well as generation and propagation of the secondaries 
is needed in order to assess in a more quantitative way observable effects of 
possible metric fluctuations on UHECRs. 

Acknowledgements We are grateful to Giuseppe Di Carlo for very useful 
and stimulating discussions. 

Note Added 

Soon after the completion of the present work, a paper [25] appeared in the preprint 
archives which addresses the problem of the fluctuations in the metric of space-time 
and the corresponding violations of LI. Our results include the results of that paper, 
but also represents a generalization to the particles with energy below the classical 
threshold for photopion production, that in our opinion imply a possibly stronger 
constrain on the effect of fluctuations of space-time on cosmic rays propagation. 
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